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Replacement of electronic switching circuits in 
computing and telecommunication systems with purely optical 
devices offers the potential for extremely high throughput 
and compact information processing systems. The potential 
application of organic materials containing molecules with 
large nonresonant nonlinear effects in this area have 
triggered intensive research during the last decade. 

Interest on this area was due to two facts: (a) that many 

organic materials show nonlinearities that are orders of 
magnitude larger than those of conventional inorganic 
materials such as lithium niobate and potassium dihydrogen 
phosphate, and (b) that organic materials much flexibility 
in terms of molecular designs. 

Some of the desirable characteristics that these 
materials should have are that they be transparent to the 
frequency of the incident laser and its second or third 
harmonic, that they have a high damage threshold and, in the 
case of second-order effects, that their crystal structure 
or molecular orientation be acentric. Since polymeric 
assemblages can enhance the nonlinear response of organic 
molecules severalfold, efforts have been directed toward the 
synthesis of thin films with interpenetrating lattices of 
electroactive molecules. The goal of this theoretical 
investigation is to predict the magnitude of the molecular 
polarizabilities of organic molecules that could be 
incorporated into films. These calculations are intended to 
become a powerful tool to assist material scientists in 
screening for the best candidates for optical applications. 

There are two possible approaches in the theoretical 
prediction of the magnitude of molecular susceptibilities of 
organic molecules: the sum-over-state approach, which is 
based on a perturbative expansion in terms of molecular 
eigenstates, and the finite-field-perturbation approach, 
which diagonalizes a molecular Hamiltonian that includes a 
static-field dipolar perturbation. In the first approach, a 
perturbation expansion in terms of excited states of the 
unperturbed molecule is used. In the second approach, 
ground state energies are calculated directly by a 
variational method in the presence of a perturbing field. 

In either case, any semiempirical or ab-initio method may be 
utilized. In general, ab-initio theories are used in 
attempting to reproduce absolute magnitudes but 
semiempirical models have been useful to detect trends among 
related molecules. For a summary of theoretical screening 
methods, please see Penn, Cardelino, Moore, Shields and 
Frazier, 1991 [5]) . 

The procedure that was developed for the present 
calculations is based on the static-field approach, and is a 
modification to the method developed by Dewar and Stewart, 
1984 [2] for calculating molecular linear polarizabilities. 
When a molecule is placed in an electric field, a given 
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component of its polarization can be expressed by the 
following summation: 

P o - n q + X «qj Fj + ? ; Pqjk Fj^lc + I Yqjkl + --- HI 

where p a represents the q component of the polarization, 
the i component of the applied electric field, H the 
permanent dipole moment, a the linear polarizability, p the 
second-order polarizability, y the third-order 
polarizability, etc. The magnitude of the vector part of p, 
corresponds to the effective measure of hyperpolarizability 
for second harmonic generation (SHG) : 

Pvec < P * 2 + Py 2 + h* > 1/2 t2! 

wJi 03 T 0 X/ y , and z tli 0 mo 1 0 cu 1 a coo3rdinat0 axss . TI 10 

symm 0 try of P with r6Sp6Ct to intercha.ncf0 of th0 last two 
indices determinos that the vector components bet 

0 q - P qqq + ( 2 P qjj + 2 I P jqj >/3 t31 

Finally, the value of the average (scalar part) third-order 
polarizability y, which is related to third-harmonic 
generation, can be expressed as: 

Y = (Yxxxx + Yyyyy + Yzzzz + 2 Y X xyy + 2 Yyyzz + 2 Yxxzz )/ 5 

Dewar and Stewart [2] method utilizes four isolated 
point charges (or sparkles, with no associated electrons or 
atomic orbitals) in the description of the molecule, 
allowing for the corresponding changes in polarization and 
gnergy to be calculated directly. The polarization and 
energy are obtained within the MNDO approximation, 
implemented in the MOPAC [4] quantum mechanical package. 

The computer coding required for the calculation of 
second-order polarizabilities was developed before the 
summer of 1991 (Cardelino et al., 1991 [1]). The original 
version of MOPAC used involved the use of fields along the 
three Cartesian coordinates. In order to obtain the 27 
required elements of the (3 tensor, calculations with fields 
along the bisector of all quadrants were added to MOPAC. In 
total, molecules are subjected to electrical fields of 14 
different strengths along each of the molecular axes and 
along the bisector line of each quadrant. Thus, MOPAC is 
used to obtain polarization data which are subsequently 
utilized by another program to perform polarization 
expansions in terms of field strengths. 

The procedure to obtain second-order polarizabilities 
was applied to nit roan i lines (Cardelino, Moore and Stickel, 
1991 [1]) and benzophenone derivatives (Moore and Cardelino, 
1991 [3]). During the summer of 1991, this code was 
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reviewed and substantially simplified. The simplified 
procedure was applied to mono-, di-, and tri-substituted 
naphthalene, quinoline and isoquinoline. A manuscript to be 
submitted for publication is under preparation. 

Most of the effort during the summer of 1991 was 
directed toward extending the procedure for second-order 
polarizabilities to obtain third-order polarizabilities. 
Based on equation [1], a polarization expansion for fields 
along one Cartesian coordinate are as follows: 

Pq = M-q + a qi F + Pqii ^ + Yqiii + . . . [5] 

where F = F^. When the fields are along the bisector of a 
quadrant, the polarization expansion can be expressed with 
two different equations, depending if the two components 
have the same or opposite signs: 

Pq = M-q + ( a qi + a qj) ^ + (Pqii + Pqjj + ^ Pcji j ^ + 

(Yqiii + Yqjjj + ^ Yqiij + ^ i j j ^ ^ *•' (^] 

where F = F^ = Fj; and, 

Pq = M-q + (®qi ~ a qj) ^ + (Pqii + Pqjj ~ ^ Pcqi j ^ ^ + 

(Yqiii ” Yqjjj “ 3 Yqiij + 3 Yqijj) ••• 

where F = F^ = -Fj. Sumation and subtraction of the third- 
order coefficients of equations [6] and [7] result in: 


2 Yqiii + 6 Yqijj 
2 Yqjjj + 6 Y qiij 


To calculate y from equation [4], the only terms that are 
needed are those in which q=i or q=j. If Kleinman symmetry 
is assumed and if q=i then qij=iij=jii and qii j=iii j= jiii; 
similarly, if q=j then qij=jij=ijj and qi j j= ji j j=i j j j . 

Then, equations [6] and [7] may be written as: 

Pi ~ M’i ~ (®ii — a ij) F - (Piii + P i j j — 3 Pjii) F 2 

(Yiiii i Yijjj ± 3 Yjiii) F 2 = + 3Yiijj F ^ [10] 

and 


Pj M- j ~ (Ciji i ®jj) F (Pjii Pjjj — 3 Pijj) F ^ 

(Yjiii — Yjjjj + 3 Yijjj) F ^ = ± 3Yiijj F^ [11] 


in which case, all terms in the left side of the equations 
can be obtained from calculations over fields along only one 
Cartesian coordinate at a time. In summary, the nonlinear 
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terms needed to obtain the iijj terms are: ijj» jii* ijjj/ 
and jiii. Figures 1 and 2 show the dependence of the 

coefficients required for the calculation of y with the ^ _ 

degree of the polynomial expansions for one particular IP 

molecule. All terms are very stable except for two 4th 

degree terms: yyyy and yxxx. These instabilities deserve 

further study. The section of the computer program that 

performs calculation of y is under way. 

Third-order polarizability calculations will be 
performed on l-pyrryl-5-x- (1, 3) -dipentyne and l-pyrryl- 6 -x- 
( i f 3 ) -dihexyne , where x will be hydroxyl, nitro, or amino 
groups. In particular, there has been a patent disclosure 
by Samuel P. McManus, Donald 0. Frazier, and Mark S. Paley 
on the hydroxyl molecule. During this summer, energy 
calculations on these molecules were also performed. 
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Fi'mrgs 1 and 2 . Values of the third and fourth degree 
coefficients of the expansions of polarization as a function 
of field strength. 
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